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ABSTRACT 


We describe our efforts in studying and comparing the ocean color data derived from 
the Japanese Ocean Color and Temperature Scanner (OCTS) and the French Polarization 
and Directionality of the Earth’s Reflectances (POLDER). OCTS and POLDER were 
both on board Japan’s Sun-synchronous A dv anced Earth Observing Satellite (ADEOS-1) 
from August 1996 to June 1997, collecting about 10 months of global ocean color data. 
This provides a unique opportunity for d evelo ping methods and strategies for the merging 
of ocean color data from multiple ocean color sensors. In this paper, we describe our ap- 
proach in developing consistent data processing algorithms for both OCTS and POLDER 
and using a common in situ data set to vica r iously calibrate the two sensors. Therefore, 
the OCTS and POLDER-measured radi ances are effectively bridged through common in 
situ measurements. With this approach in processing data from two different sensors, the 
only differences in the derived products from OCTS and POLDER are the differences 
inherited from the instrument characteri stics. Results show that there are no obvious bias 
differences between the OCTS and PO LDER -derived ocean color products, whereas the 
differences due to noise, which stem fro m v ariations in sensor characteristics, are difficult 
to correct. It is possible, however, to re duce noise differences with some data averaging 
schemes. The ocean color data from OC T S a nd POLDER can therefore be compared and 
merged in the sense that there is no significant bias between two. 


1. Introduction 


The Japanese Ocean Color and Temperature Scanner (OCTS)i and French Polariza- 
tion and Directionality of the Earth’s R eflec tances (POLDER)^ were both flown on the 
Japanese Sun-synchronous Advanced Eartb Observing Satellite (ADEOS-1) from August 
1996 to June 1997, collecting about 10 months of global ocean color data. ADEOS-1 
was on a polar orbit at an altitude of 800 km with local crossing time (descending node) 
at around 10:40 AM. This was the first time in history that two ocean color sensors were 
on board the same platform and viewe d th e global ocean with the same temporal and 
similar global spatial coverages. Therefore, it provides an ideal case for studying and 
comparing ocean color data derived fro m two different sensors, hence allowing develop- 
ment of a strategy for ocean color data merger from multiple ocean color sensors. The 
primary goal of the NASA Sensor Intercpmparison and Merger for Biological and Inter- 
disciplinary Oceanic Studies (SIMBIOS) pr oject is to develop methods for meaningful 
comparison and merging of data products from multiple ocean color missions.^ In a re- 
cent study, Wang and Franz^ show that, using a vicarious inter-calibration approach be- 
tween the modular optoelectronic scanner (MOS)^ and the Sea-Viewing Wide Field-of- 
View Sensor (SeaWiFS)^'^ (SeaWiFS ocean color data used as “truth”), the ocean optical 
property data derived from two sensors eEn be meaningfully compared. The bias differ- 
ences are significantly reduced between t he t wo products. In that study measurements 
from the two sensors had a temporal difference of about 90 minutes. In this paper, we 
compare the ocean color data derived fr om OCTS and POLDER measurements using 
consistent data processing algorithms for both sensors and vicarious calibrations based on 
a common in situ data set from the Marine O ptical Bouy (MOBY).® Therefore, differ- 
ences in the derived ocean color products from OCTS and POLDER are primarily associ- 
ated with differences in instrument char acteris tics. We first give a brief overview of the 
sensor characteristics of the OCTS and the P OLDER and their differences. Detailed al- 


1 


2 


gorithm descriptions for processing OCXS and POLDER data are then presented. Next, 
we outline the vicarious calibration scheme in which the OCXS and POLDER measured 
radiances are calibrated using a common in situ MOBY data set. Finally, we provide re- 
sults comparing OCXS and POLDER oceancolor data with global in situ measurements, 
and compare a series of OCXS measure ments with those of POLDER over the Sargasso 
Sea and the Bermuda area. 

2. Instrument Characterizations 

OCXS, which was built and operated by the National Space Development Agency 
(NASD A) of Japan, is an optical radio meter ^ with a scanning swath of 1400 km and a 
spatial resolution at nadir of 0.7 km. O CX S has 8 spectral bands (out of 12) for ocean 
color measurement centered at wavelengths 412, 443, 490, 520, 565, 670, 765, and 865 
nm. Xhe bandwidths of the first 6 bands are 20 nm and those of the last two near-infrared 
(NIR) bands are 40 nm. Xhe NIR bands can be used for the atmospheric corrections.^'^^ 
In addition, OCXS has a capability of operationally tilting the sensor ±20° away from na- 
dir to minimize sun glint contamination. In c omparison. POLDER, which was built and 
operated by the Centre National d’Etudej^patiales (CNES) of France, is a multi-angle 
imaging radiometer/polarimeter with a 242x274 CCD pushbroom matrix array. ^ 
POLDER has a wide field of view (2400 1^ a spatial resolution at nadir of ~6-7 km, 
and 6 spectral bands centered at 443, 490, 565, 670, 765, and 865 nm for ocean color 
measurement. Xhe bandwidths of the first fou r visible bands are 20 nm and those of the 
two NIR bands are 40 nm. Both 670 and 865 nm bands are designed to measure the po- 
larized radiance in addition to another 4 4Xirm polarized band. For a given earth ground 
scene (target), POLDER can provide up to 14 near-contemporaneous measurements at 
various viewing directions, therefore it is u s eful for deriving atmosphere and/or ground 
bi-directional reflectance distribution f u ncti o ns (BRDF). Obviously, these two sensors 
have many different design characteristics,“even though their main goal of retrieving the 
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global ocean color products is the same. For that purpose, they have 6 common spectral 
bands at 443, 490, 565, 670, 765, and 86 5 nm . 

3. Algorithm Descriptions 

The data processing algorithms used for OCXS and POLDER are described in the 
following two sections. The processing algorithms are implemented in the Multi-Sensor 
Level- 1 to Level-2 (MSL12) software whic h is freely distributed as part of SeaWiFS 
Data Analysis System (SeaDAS) software package.'^ 

A. Data Processed to Level- IB 

The OCXS and POLDER Level- IB data are usually obtained from Japan’s NASDA 
and the French ONES data centers. Ho weve r, during the ADEOS-I mission lifetime, ap- 
proximately 450 GB of real-time OCXS data were collected by the SeaWiFS project 
through the National Oceanic and Atmospheric Administration (NOAA) ground stations 
at Wallops, Virginia and Fairbanks, Alaska. Since standard OCXS Level-0 to Level-IB 
software was not publicly available, the SIMBIOS project engaged in developing and en- 
hancing the Level-0 to Level- IB processing software, with specific attention given to im- 
proving geolocation and image registration. Herein, we briefly describe the OCXS Level- 
0 to Level- IB processing software developed by the SIMBIOS project. 

In general, the processing of Level-0 data to Level-IB includes conversion of raw 
sensor counts to physical units and the assignment of geolocation information to each ob- 
servation (pixel, detector, and band). Several factors in the design of the OCXS instru- 
ment complicate this process. Each OC XS band is divided into 10 individual detectors, 
where each detector is associated with a separate scan line. Each detector has a slightly 
different and apparently non-linear responSivity, which gives rise to horizontal striping in 
the OCXS imagery. 13 The basis for the SIMBIOS calibration is the NASDA Version 3 
preflight calibration and relative detector calibration, i'* The relative detector calibration 
substantially reduces the striping effects, but low-level artifacts are still visible and they 


4 


can become significant following removal of the large, comparatively smooth atmos- 
pheric signal. 

The design of the OCTS scan mirr o r an d focal plane results in severe misalignment 
of the line-of-sight for different spectral band s viewing at essentially the same time.i^ To 
perform atmospheric correction or derive chlorophyll concentration at a given location 
requires that the bands be co-registered to that location. The SIMBIOS project developed 
a registration technique in which an idealized focal-plane geometry is defined, and all 
bands are resampled using a nearest neighbor approach to match the idealized ground 
track. This method is similar in concept to the approach of using OCTS band 1 as a near- 
est neighbor matching target, 13 but it also eliminates scan to scan overlap. Furthermore, 
while the registration process can correct for differences in geolocation between bands, 
the resampled observations come from different lines-of-sight with unique sensor-to- 
ground and Sun-to-ground geometries. To ensure accurate atmospheric correction, the 
true solar and sensor view angles for eac h ban d are saved in the Level- IB file for use in 
the Level- IB to Level-2 processing. Note that this is not possible with NASD A standard 
Level- IB products because only nominal angles are saved. 

Before the registration process c an begin , however, the observations from each band 
must be accurately navigated. The geolocation algorithm used for SIMBIOS OCTS proc- 
essing is a modified version of the meth od used for SeaWiFS navigation.15 tech- 

nique has the potential to yield an exact solution to the geolocation of each scan line 
when the spacecraft position and sensor geometry are accurately known. To improve 
geolocation accuracy, the SIMBIOS project derived time-dependent attitude and tilt ad- 
justments using an automated navigation assessment technique based on island targets. 
Results show that, using the SIMBIOS Level- IB data processing, the OCTS navigation is 
accurate to approximately one kilometer. 

The OCTS Level- IB data processed with the SIMBIOS code usually agree well with 
the NASDA standard Level- IB product.~Fo'r OCTS bands 1-8, the radiance frequency 
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distributions from the two processes are nearly identical. However, the radiance differ- 
ences based on a pixel by pixel comparisori can disagree due to small differences in navi- 
gation and band coregistration, i.e., in such comparisons, the corresponding pixels from 
the NASD A standard Level- IB and the SIM BIOS Level-0 to Level- IB routine are usu- 
ally not the same. In a typical example, these differences were found to be 0.54%, 
0.53%, 0.54%, 0.85%, 1.00%, 1.62%, 2.20%, and 2.67% for bands 1-8, respectively. 

B. Data Processed from Level- IB to Levgl-2 

For ocean color remote sensing. Level- IB to Level-2 data processing mainly in- 
volves the atmospheric correction,9‘ii4^48 in which the radiance contributions from the 
atmosphere and ocean surface as well as their interactions are removed from the sensor- 
measured signals at the top of the atmosphere (TOA). The radiance contributions from 
the ocean, i.e., the water-leaving radiances just above the sea surface, can then be re- 
trieved. A simple empirical algorithm^^' ^^ can then be used to derive the ocean chloro- 
phyll concentration from the normalize d wa ter-leaving radiances. In the following, a 
consistent atmospheric correction algorithm applied to OCXS and POLDER as well as a 
bio-optical algorithm tuned to OCXS and POLDER spectral bands are briefly described. 
The data processing procedures, in particular for POLDER, are outlined. 

7. Atmospheric Corrections 

Based on a study by Wang,22 the SIMBIOS project has developed the MSL12 soft- 
ware for processing data from multiple ocean color sensors. MSL12 is an implementa- 
tion of the standard SeaWiFS atmospheric correction algorithm,^47 vvith some generali- 
zations to allow for application to OCXS, POLDER, and other ocean color sensors (e.g., 
MOS).^'22 Therefore, a consistent atmospheric correction algorithm can be used to proc- 
ess data from OCXS and POLDER. 

Similar to SeaWiFS, both OCXS and POLDER have two NIR bands at 765 and 865 
nm for the atmospheric correction. However, the 765 nm spectral band completely en- 
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compasses the oxygen A-band absorption. Figure 1 compares the OCTS and POLDER 
765 nm band spectral response functions (SRF) with that of SeaWiFS. Clearly, both 
OCTS and POLDER 765 nm band SRFs are almost identical to that of SeaWiFS. There- 
fore, the same SeaWiFS oxygen A-band absorption correction algorithm's can be used 
for both OCTS and POLDER. A small difference in the SRFs (mainly between 
POLDER and SeaWiFS) is within the un certa inty of the correction algorithm. ^-24 

There are some major improvements in MSL12 for the OCTS and POLDER data 
processing which are concurrent with changes in processing algorithms for the SeaWiFS 
3rd data reprocessing in May 2000.^5 Some specific atmospheric correction algorithm 
improvements include: 

1 . An iterative correction algorithm, which accounts for the ocean radiance contri- 
butions at two NIR wavelengths (765 and 865 nm), was developed and implemented in 

the data processing. ^6 

2. Improved Rayleigh radiance lookup table, in which the Rayleigh (air molecules) 
scattering radiance is a function of the sea su rface wind speed (in addition to solar and 
sensor viewing geometry and wavelength dependence),^^ was generated and imple- 
mented in the MSL12 data processing system^® for OCTS and POLDER. 

3. A sun glint contamination corr ectio n algorithm, which effectively removes the 
sun glint residual contamination radiance around the sub-solar point, was developed and 
implemented in MSL 12.29 

4. Correction of the ocean whitecapTadiance contributions^ 2 ^vas updated^S using 

some current in situ measurements. 30-32 

5. The current 12 aerosol models used in generating aerosol lookup tables for the 
atmospheric correction are: Oceanic wi th rela tive humidity (RH) of 99%, Maritime with 
RH of 50%, 70%, 90%, and 99%, Coastal with RH of 50%, 70%, 90%, and 99%, and 
Tropospheric with RH of 50%, 90%, and 99%. 25 The Oceanic, Maritime, and Tro- 


pospheric models are from Shettle and Fenn,33 whereas the Coastal model was derived 
using data from the Shettle and Fenn mo dels. ^ 

6. A correction algorithm for the sensor spectral bandpass effects on the derived 
normalized water-leaving radiances for SeaWiFS^^ was developed. Similar algorithms 
for both OCXS and POLDER were also iniplemented in MSL12. Similar to Wang et. al. 
(Figure 3 (b)),34 Figures 2(a) and 2(b) provide the spectral bandpass correction factor 
r(A) as a function of the sensor-derived (uncorrected) ratio values between the normal- 
ized water-leaving radiances at wavelengths 490 and 565 nm. Note that r(A) >1 indi- 
cates under-estimation of the derived no rmalized water-leaving radiance [X^(A)]^, 
whereas r(A) < 1 indicates over-estimationZ Fig- 2(a) shows the correction factors for the 
OCTS’s spectral bands, while Fig. 2(b) shows those for POLDER. Fig. 2(a) shows that, 
without corrections, the OCTS-derived [L^(A)]^ at 565 nm is under-estimated (~8%), 
while the values at 520 nm are over- estima ted (up to -17%) for clear ocean waters. 
Other OCXS bands (412, 443, and 490 nm) have almost no spectral bandpass effects. On 
the other hand, without corrections, the POLDER 490 nm band has ~5% under- 
estimation for the derived value in clear ocea n waters, while for the 565 nm band there is 
under-estimation from about 3.5%-1.5% o f the derived values for ocean waters ranging 
from productive to non-productive (clear) cases. Clearly, corrections of the sensor spec- 
tral bandpass effects are necessary in order to have a meaningful Level-2 product com- 
parison derived from two different sensors, e g., OCXS and POLDER. 

2. Bio-Optical Algorithms 

In this study, we use the bio-optical algorithm of ocean chlorophyll 2 (OC2)2i for de- 
riving the ocean chlorophyll concentration from the sensor-retrieved normalized water- 
leaving radiance [E„,(A)]^. Xherefore, a consistent bio-optical algorithm is used in de- 
riving chlorophyll concentration for both OCXS and POLDER. Xhe OC2 algorithm was 


modified by O’Reilly and Maritorena (per. comm.) to account for the effects of the spec- 
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tral band differences between the OCTS and POLDER 565 nm band and the SeaWiFS 
555 nm band, i.e., 


- jq(o.3164-2.I32/?+0,6303/?Vo.004«^) 


0.0708, 


( 1 ) 


where I 

fi = [p.{«0)y[p.(565)]„ (2) 

is the ratio between the normalized water - l e aving reflectances at 490 and 565 nm. The 
normalized water-leaving reflectance, is defined as 

[p.(A)],=<r[i.(A)]„/F,(A), (3) 

where [L^(A)]^ is the sensor-derived normalized water-leaving radiance and Fq(A) is the 
extraterrestrial solar irradiance. 


3. A Specific Data Processing Procedure for POLDER 

As we discussed in previous sections, for a given location (scene), POLDER can ac- 
quire as many as 14 near-contemporaneous measurements. Therefore, up to 14 values of 
the Level 2 product for the same ocean scene c an be derived. However, it was noted that, 
for a given scene, the POLDER derived_products, e.g., normalized water-leaving radi- 
ances, often have high variations (up to an order of magnitude) between various 
POLDER viewing directions. These var iatio ns are possibly related to the POLDER sub- 
pixel cloud contamination and/or band no n- linear responses.^^ It was also found that, for 
a given scene, there was no obvious cor re la t ion between the POLDER-derived 
and its viewing direction (angle). Therefore, we used an averaging scheme to eliminate 
points with high variation and to reduce data noise in the derived product. First, for a 
given pixel, the derived normalized water-leaving radiances from all possible POLDER 
viewing directions were averaged. Data points which were within 50% variation of the 
first averaged value (i.e., 0.5- 1.5 of average value) were then used to derive the final 
mean value and considered as the POLD ER -derived normalized water-leaving radiances 
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for the pixel. This averaging scheme was used in deriving the POLDER vicarious gain 
coefficients, as well as the POLDER Level 2 products. 

4. Vicarious Calibrations with MOB Y in situ Measurements 

It is well known that, for ocean color remote sensing, an on-orbit vicarious calibra- 
tion of sensor and algorithms is necessary.3^ Vicarious calibration techniques have been 
successfully applied to SeaWiFS^^'^® and also used to inter-calibrate MOS and SeaWiFS 
data.^ In this study, a vicarious calibration for OCTS and POLDER using the in situ 
measurements from the MOBY® was carried out. The sensor vicarious gain coefficients, 
A), are defined through 

( 4 ) 

where is sensor-measured Level- IB radiance from OCTS and POLDER and 

is vicariously calibrated radiance which is used in the Level-1 B to Level-2 algo- 
rithms as the correct value of L,(A). The vicarious calibration procedure was proceeded 
as follows. First, for the 865-nm band, it was assumed that gain coefficients were not 
changed. They were equivalent to the values suggested by NASDA and CNES. The in- 
flight gains for OCTS and POLDER 865-nm band were derived by NASDA and CNES 
as 0.89 and 1.05, respectively. Next, the gain coefficient of the 765 nm band was de- 
rived by forcing the atmospheric correction algorithm to select a maritime aerosol type 
over the MOBY site. The procedure was similar to the methods used in adjusting the 
SeaWiFS 765 nm band.®^ Finally, with the derived 765 nm gain coefficient, the MOBY 
in situ measurements were used to calibrate the visible bands,®® i.e., the OCTS and 
POLDER visible bands were calibrated such that, with the atmospheric corrections, the 
derived normalized water-leaving radiances were equivalent to the MOBY in situ meas- 


urements. 
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The location of the MOBY site is usually in stable dear-ocean waters and the aero- 
sols are predominantly marine. The MOBY program has been providing consistently 
high quality dear-ocean optical data since 1997. Unfortunately, there were some prob- 
lems with the early deployment of MOBY in 1996, limiting the number of usable cali- 
bration samples to five measurements spanning the period of November 1996 to February 
1997. Specifically, five in situ measurements taken on November 26 of 1996, January 
13, January 17, February 16, and February 23 of 1997 were used. Table I provides the 
derived gain coefficients for the OCTS an d P OLDER spectral bands. 

5. Results 

Using the derived vicarious gain coefficients, OCTS and POLDER data were proc- 
essed to Level-2 products from either OCTS Level-0 (data acquired from Wallops, Vir- 
ginia) and/or OCTS Level-IB data and POLDER Level-IB data using MSL12. As dis- 
cussed in the previous sections, the processing algorithms for both OCTS and POLDER 
are essentially the same as for SeaWiFS although some specific modifications are neces- 
sary depending on sensor characteristics, In the following four sections, we provide 

OCTS and POLDER matchup analyses comparing satellite retrievals with in situ meas- 
urements, ocean color data comparisons between the two sensors, as well as some data 
analyses and discussions. 

A. Global in situ Matchup Analyses 

An extensive set of in situ data take n during the time of the OCTS and POLDER 
missions was archived in the SeaWiFS Bio-optical Archive and Storage System 
(SeaBASS) data base. These data were used for OCTS and POLDER matchup analyses. 
In these comparison analyses, the satellite and in situ data were screened and obtained 
through the following procedure: (i) the differences of time and solar zenith angle be- 
tween satellite and in situ measurement s w ere less than 4 hours and 15°, respectively; 
(ii) the solar zenith angle and sensor viewing angle for OCTS and POLDER were < 60°; 
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(iii) the in situ measurement was between 9 AM to 3 PM local time; (iv) in the matchup 
analyses, the OCTS [^h-(^)]jv derived by averaging data in a 5 km circle centered at 
the location of the in situ measurement, and only data for which a coefficient of determi- 
nation (the ratio of standard deviation to average) was less than 0.5 were considered; and 
(v) as discussed previously, POLDER results were derived by first averaging over all 
possible viewing directions (up to 14 values) and data were then re-averaged after elimi- 
nating data points whose differences compared with the first averaged value were > 
±50%. 

Figures 3 and 4 provide matchup results for the OCTS and POLDER compared with 
the in situ measurements. Figs. 3(a) and 3(b) show matchup results for the normalized 
water-leaving radiance derived from OCTS NASDA standard Level- IB data 

and Wallops Level-0 data (using SIMBTOS Level-0 to Level-IB processing), respec- 
tively. Fig. 4(a) shows the matchup results derived from POLDER CNES Level- IB data 
compared to the in situ measurements. Figs. 3(c) and 3(d) compare the OCTS derived 
two-band ratio value [L^(A)]^/[L^,(565)]^ from NASDA Level-IB and Wallops Level-0 
data, respectively, with results from in sit u m easurements, while Fig. 4(c) shows the same 
results from POLDER matchup comparisons. The matchup analyses for the two-band 
ratio value in the derived are useful since they are inputs for the ocean color 

bio-optical algorithms. 

Results in Figs. 3 and 4 show that, after vicarious calibration using the MOBY data, 
the derived [Av('^)];y,/[^B-(565)]^ from the OCTS and POLDER measure- 

ments compared reasonably well with the in situ data. Note that each plot gives the linear 
fit showing coefficients of slope, intercept (Int), and the correlation coefficient (R). (A 
perfect match would have slope of 1, intercept of 0, and the correlation coefficient of 1). 
It appears, however, that the OCTS-derived values were slightly over-estimated 

(slope > 1.0), while the POLDER results were slightly under-estimated (slope < 1.0). 
The derived two-band ratio values [l^(A)]^/[l,^(565)]^ compared better to the in situ 
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measurements than those of [l„(A)]^, indicating that the derived ocean chlorophyll con- 
centrations are reasonably accurate. Results also indicate that the POLDER-derived 
products are more noisy than the OCXS data since the correlation coefficients derived 
from the POLDER matches are lower than those from the OCXS comparison. 

B. OCXS and POLDER Comparisons 

A series of measurements from O CXS and POLDER over the Sargasso Sea and the 
Bermuda area were collected for comparison. Xhe satellite observations were obtained 
from a co-located 3°x3° box over the Sargasso Sea (latitude of 23.5°N to 26.5°N and 
longitude of 70°W to 73°W) and the Bermuda area (latitude of 30°N to 33°N and longi- 
tude of 62.7°W to 65.7°W). Over the Sargasso Sea and the Bermuda area, OCXS data 
were all collected through the ground stati on at Wallops, Virginia. Xherefore, OCXS data 
were processed from Level-0 to Level- IB u sing the SIMBIOS code, and then to Level-2 
using MSL12. Figure 5 provides color images for the OCXS and POLDER derived 
comparisons for wavelengths of 443, 490, and 565 nm over the Sargasso Sea 
and the Bermuda area. Fig. 5(a) is for the results over the Sargasso Sea, while Fig. 5(b) is 
for comparisons over the Bermuda area. Xhe color scales for images are indi- 

cated at the bottom of Fig. 5(b). [L^(A)|^ for 443 and 490 nm is scaled from 0 to 3 (mW 
cm'^ |im’’ sr'*), while for 565 nm is scaled from 0 to 0.8 (mW cm'^ |im'' sr"'). 

In Fig. 5, columns 1-3 are the OCXS-derived images for 443, 490, and 565 nm, 

and columns 4-6 are the POLDER-derived results for 443, 490, and 565 nm. In addition, 
at the left side of the images, there are indicated dates when OCXS and POLDER meas- 
urements were acquired. Over the Sargasso Sea, OCXS and POLDER data were acquired 
on November 20 and December 19 of 1996, and January 26, February 17, March 11, 
April 14, and May 20 of 1997 (there were data from each month for seven months), while 
data were acquired on December 12 of 1996, and March 19, May 3, May 7, and June 10 
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of 1997 over the Bermuda area. In Fig. 5, each image represents a 3°x3° box over the 
Sargasso Sea or the Bermuda area. 

Results in Fig. 5 show a reasonable agreement between OCTS and POLDER derived 
[L„(A)]^ values. Over the Sargasso Sea, both spatial and temporal variations of ocean 
color were well captured by both OCTS and POLDER. For example, [L„,(A)]^ values at 
443 nm changed from November (note the spatial variations in the November image from 
both OCTS and POLDER) to a lower value in December and January, then elevated in 
February, and finally reached their highest values in the spring. The spatial [T„,(A)]^ 
distributions from OCTS and POLDER are comparable. However, there are some obvi- 

ous differences, e.g., POLDER has a larger pixel size (low spatial resolution) and it ap- 
pears that the POLDER results show more noisy than results from OCTS. Some spar- 
kling points in the POLDER [T„(A)]^ images are mostly cloud-edge effects, probably 
due to sub-pixel cloud contamination because of POLDER low spatial resolution meas- 
urement. 

Figure 6 compares the histogram (%) of the retrieved [L„(A)]^ between OCTS and 
POLDER over the Sargasso Sea and the Bermuda area. Figs. 6(a)-6(g) are comparison 
results for the Sargasso Sea for data acquired from November 20, 1996 to May 20, 1997, 
while Figs. 6(h)-6(l) are for the Bermuda area for measurements acquired from December 
12, 1996 to June 10, 1997. The OCTS and POLDER [L,,(A)]^ at 490 and 565 nm usu- 
ally compare quite well, with small variations between the two measurements, while the 
POLDER [L„,(A)]^ at 443 nm shows significant variations compared with the OCTS re- 
sults (sometimes high and sometimes low). The shape of the histogram distributions 
from OCTS and POLDER, however, are generally similar. Table 2 gives quantitative 
comparisons of the retrieved [L„(A)]^ between OCTS and POLDER for the mode (peak) 
values in the histograms, as presented in Fig. 6. Table 2(a) and 2(b) are results for the 
Sargasso Sea and the Bermuda area, respectively. Although there are some variations in 
the derived [T„,(A)]^ values between OC TS and POLDER, the average ratio values of 
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OCTS/POLDER for 443, 490, and 565 nm for all cases in Table 2 are: 1.0005, 

1.012, and 0.9831. These results indica te that, with our approach in processing OCTS 
and POLDER data, there are no obvious bias differences in the derived [L„,(A)]^ between 
the two sensors, whereas the noise differences are difficult to correct since they depend 
mainly on the sensor characteristics. 

C. OCTS and POLDER Chlorophyll-a Comparisons 

Figure 7 provides comparisons of chlorophyll-a concentration derived from OCTS 
and POLDER using the OC2 algorithm, which was modified for OCTS and POLDER, 
over the Sargasso Sea and the Bermuda apa. Fig. 7(a) shows results for the Sargasso 
sea, while Fig. 7(b) shows those of the Bermuda area. The color images are scaled loga- 
rithmically from 0.01 to 5 (mg m'^), and the color scale is indicated at the bottom of Fig. 
7(b). In Fig. 7, columns 1 and 2 are color images, respectively, for the OCTS and 
POLDER-derived chlorophyll-a concentrations, while column 3 shows a histogram com- 
parison of the distributions between the two color images. 

Over the Bermuda area (Fig. 7(b)), there were some in situ measurements which 
were acquired at about the same time as for some the OCTS and POLDER measure- 
ments. These data are archived in the Se aB ASS data base. The in situ data were taken 
near the center of the Bermuda 3°x3° b ox are a (center of the images). These in situ data 
are shown as vertical bars in the histograms in Fig. 7(b). Therefore, for the case of De- 
cember 12, 1996, the OCTS and POLDER-derived chlorophyll values appear to be un- 
der-estimated, as compared with the in situ data. However, there were no retrievals from 
OCTS and POLDER over the center and top-most two-thirds part of the image box area. 
Therefore, the values from the satellites may not be representative for comparison with 
the in situ measurements. The same is als o true for the case of June 10, 1997, where 
there were no retrievals for the bottom- most two-thirds part of the images. On the other 
hand, when coverage did coincide with field data as for both March 19,1 997 and May 7, 
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1997 both satellites compared well with in situ measurements. Note that in situ data for 
May 7, 1997 and June 10, 1997 were obtained one day earlier (May 6, 1997) and one day 
later (June 11, 1997), respectively. 

There are obvious variations in the derived chlorophyll-a concentrations between 
OCXS and POLDER. However, it is fou n d that, for the chlorophyll-a concentration 
comparison, the average ratio value (mode in histogram) of OCTS/POLDER for all the 
Sargasso Sea and the Bermuda area cases (total of 12) is 0.9907. 

D. Discussion 

In examining the results of the comparisons between OCXS and POLDER, it is clear 
that there are no bias differences in the derived products between the two sensors. How- 
ever, it has been consistently shown that, compared to OCXS, the POLDER-derived 
ocean color data are often more noisy, in particular, for data derived from the POLDER 
443-nm band. One possible reason for such variability could be due to POLDER sensor 
non-linearity .35 Results have shown35 t hat t he POLDER sensor non-linear responses un- 
der-estimate radiance values at the loweLradiance ranges, while they over-estimate val- 
ues for the high radiance ranges. It was reported35 that the POLDER 443 nm band has 
about a 3% non-linear variation from low to high radiances for ocean color measure- 
ments. Other bands also have non-linear responses, but with much smaller magnitudes.35 

Some high sparkling values in the POLDER-derived ocean color products, however, 
are mostly from the cloud-edge effects, due to sub-pixel cloud contamination. With a 
low spatial resolution pixel, POLDER measurements are more likely to be contaminated 
by sub-pixel clouds. The cloud contaminated pixels could be mis-identified as clear be- 
cause their radiance values are not high enough to trigger the cloud threshold. 

6. Conclusion 

We describe a procedure to vicariously calibrate OCXS and POLDER Level- IB data 
using the MOBY in situ measurements. It effectively bridges the OCXS and POLDER 
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Level- IB data through the use of a common set of MOBY in situ measurements. The 
vicarious calibration assumes the band M5 nm gain coefficients are not changed, while 
the visible bands are adjusted such that t he derived normalized water-leaving radiances 
are forced to be equivalent to the MOBY in situ measurements. Using the derived gain 
coefficients, a consistent atmospheric correction algorithm is used to process OCXS and 
POLDER data from Level- IB to Level-2 products. Therefore, the differences in the 

OCXS and POLDER-derived ocean color products are mainly from the differences of the 

sensor characteristics, i.e., with this approach, the retrieved products would be identical if 

the two sensors were identical. 

We demonstrate the efficacy of vicariously re-calibrating OCXS and POLDER with 
a common in situ data set and using a consistent atmospheric correction algorithm for 
data processing. Both the OCXS and POLDER-derived normalized water-leaving radi- 
ances compare reasonably well with independent in situ data. Over the Sargasso Sea and 
the Bermuda area, the ocean color products derived from OCXS have a good agreement 
with data fi-om POLDER measurements. It is particularly important that there are no bias 
differences between the OCXS and PO LDE R-derived ocean color data, i.e., the overall 
ratio values between the OCXS and P OLD ER-derived normalized water-leaving radi- 
ances at 443, 490, and 565 nm are: 1.0 00 5 , 1.012, and 0.9831, while it is 0.9907 for the 
ratio value in the derived chlorophyll-a concentration. However, the noise differences, 
which are usually inherited from sensor characteristics, are difficult to remove. These 
noise differences can be possibly reduced using data averaging schemes (e.g. averaging 
Level-2 data to Level-3 products). It is found that the POLDER-derived [L„,(A)]^ are 
usually more noisy than those of OCXS, in particular, for the 443 nm band products. 
This most likely resulted from the POLDER sensor non-linearity response problem.35 On 
the other hand, effects of the sub-pixel c loud contamination around cloud edge lead to 
high sparkling values in the POLDER-retrieved ocean color products. 
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The OCXS and POLDER data processing procedures and algorithms are planned to 
be implemented in MSL12 within the SeaDAS software , 12 which is freely available. 
Therefore, interested scientists can process OCXS and POLDER data using the widely 
distributed SeaDAS software package. 
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Figure Captions 

Figure 1. The OCXS and POLDER 765 ran band spectral response functions compared 
with that of SeaWiFS. 

Figure 2. Spectral bandpass correction fact or as a function of the sensor-derived 
two-band ratio values in the normalized w ater-leaving radiances between wavelengths 
490 and 565 nm for (a) the OCXS bands T12, 443, 490, 520, and 565 nm and (b) the 
POLDER bands 443, 490, and 565 nm. The curves are the lest-squares fit. 

Figure 3. The OCXS-derived ocean color parameters compared with the in situ meas- 
urements for (a) [Lj„(A)]^ with OCXS NASDA Level-IB data; (b) with OCXS 

Wallops Level-0 data; (c) ratio [L„(A)]^/[ 1^(565)]^ with OCXS NASDA Level-IB data; 
and (d) ratio [L,,(A)]^/[L,^(565)]^ with OCXS Wallops Level-0 data. Each plot indicates 
a linear fit coefficients of slope, intercept (Int), and the correlation coefficient (R). 

Figure 4. The POLDER-derived ocean color parameters compared with the in situ 
measurements from POLDER CNES Lev el-IB data for (a) [L^(A)]^ and (b) ratio 
[L^(A)]^/[L,^(565)]^ values. Each plot indicates a linear fit coefficients of slope, inter- 
cept (Int), and the correlation coefficient (R). 

Figure 5. The OCXS-retrieved norma liz ed water-leaving radiance images at wave- 
lengths 443, 490, and 565 nm compared with those derived from POLDER measurements 
for ocean water over (a) the Sargasso Sea from Nov. 20, 1996 to May 20, 1997 and (b) 
the Bermuda area from Dec. 12, 1996 to June 10, 1997. 

Figure 6. The histogram (%) of the OCXS-retrieved normalized water-leaving radiances 
at wavelengths 443, 490, and 565 nm compared with those from POLDER measure- 
ments as in Fig. 5 for (a)-(g) the Sargasso Sea from Nov. 20, 1996 to May 20, 1997 and 
(h)-(i) the Bermuda area from Dec. 12, 1996^fo June 10, 1997. 
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Figure 7. The OCTS-retrieved chlorophyll-a concentration compared with those derived 
from POLDER measurements for (a) the Sargasso Sea from Nov. 20, 1996 to April 14, 
1997 and (b) the Bermuda area from Dec. 12, 1996 to June 10, 1997. Some of the in situ 
measurements are also indicated as vertical bars in the histograms over the Bermuda area 
in(b). 
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Table 1 . The derived vicarious gain coefficients for OCXS and POLDER spectral 
bands. 


Wavelength 

Gain Coefficient 

(nm) 

OCXS 

POLDER 

412 

1.12426 

— 

443 

1.01539 

1.06465 

490 

0.95084 

1.02350 

520 

1.01784 

— 

565 

1.03255 

0.97541 

670 

1.00859 

1.01845 

765 

0,92093 

1.02946 

865 

1.00000 

1.00000 
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Table 2(a). The OCTS/POLDER comparisons in the retrieved (mode 

value) over the Sargasso Sea from Nov. 20, 1996 to May 20, 1997. 


Date 



[^■-(565)]/ 

OCTS 

POLDER %Diff 

OCTS 

POLDER %Diff 

OCTS 

POLDER %Diff 

Nov. 20, 1996 

1.8015 

1.7500 

2.90 

1.2209 

1.1503 

5.95 

0.1975 

0.2844 

-36.07 

Dec. 19, 1996 

1.5489 

1.4115 

9.28 

1.0491 

1.1101 

-5.65 

0.1519 

0.1606 

-5.54 

Jan. 26, 1 997 

1.6381 

1.7798 

-8.29 

1.1048 

1.0874 

1.60 

0.1643 

0.1587 

3.46 

Feb. 17, 1997 

1.8567 

1.7216 

7.55 

1.1591 

1.1659 

-0.58 

0.1841 

0.2048 

-10.62 

Mar. 11, 1997 

1.8705 

2.2018 

-16.27 

1.1722 

1.2684 

-7.89 

0.1717 

0.2077 

-18.98 

Apr. 14, 1997 

2.0278 

2.0893 

-2.99 

1.2324 

1.2880 

-4.41 

0.1872 

0.1877 

-0.22 

May 20, 1997 

1.9583 

2.0156 

-2.88 

1.1064 

1.1901 

-7.29 

0.1576 

0.1525 

3.29 


tunit of (mW cm'^ pm'' sr’’). 


Table 2(b). The OCTS/POLDER comparisons in the retrieved [L^(A)]^ (mode 
value) over the Bermuda area from Dec, 12, 1996 to June 10, 1997. 


Date 


[i.(490)]„t 



OCTS 

POLDER 

%Diff 

OCTS 

POLDER %Diff 

OCTS 

POLDER %Diff 

Dec. 12, 1996 

1.4346 

0.9627 

39.37 

1.0865 

0.9672 

11.62 

0.2093 

0.1777 

16.33 

Mar. 19, 1997 

1.4925 

1.5656 

-4.78 

1.0410 

1.0329 

0.78 

0.2104 

0.2530 

-18.33 

May 3, 1997 

1.6149 

1.9393 

-18.25 

1.1171 

1.1603 

-3.79 

0.1616 

0.1515 

6.45 

May 7, 1997 

1.7566 

1.8485 

-5.10 

1.1937 

1.0473 

13.07 

0.2121 

0.1515 

33.33 

June 10, 1997 

1.3194 

1.5309 

-14.84 

1.1024 

1.0383 

5.99 

0.1616 

0.1818 

-11.76 


tunit of (mW cm'^ pm"' sr'*). 
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SIGNIFICANT RESULTS 

We describe our efforts in studying and comparing the ocean color data derived from 
the Japanese Ocean Color and Temperature Scanner (OCTS) and the French Polarization 
and Directionality of the Earth’s Reflectances (POLDER). We describe our approach in 
developing consistent data processing algorithms for both OCTS and POLDER and using 
a common in situ data set to vicariously calibrate the two sensors. Therefore, the OCTS 
and POLDER-measured radiances are effectively bridged through common in situ meas- 
urements. The ocean color data from OCTS and POLDER can therefore be compared 
and merged in the sense that there is no significant bias between two. 
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